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Introduction  and  summary . Maximum  likelihood  estimates  of  the 
parameters  of  a bivariate  normal  distribution  are  obtained  for  a sample 


in  which  only  those  observations  falling  in  a specific  region  can  be 
measured,  all  other  observations  being  called  "unmeasured  observations". 
Two  cases  are  treated,  the  number  of  unmeasured  observations  being  un- 
known (Case  I)  or  known  (Case  II),  Explicit  expressions  are  obtained 
when  the  region  of  truncation  is  a rectangle  or  an  infinite  strip.  The 
asymptotic  covariance  matrix  is  obtained  simultaneously  with  the  solution. 

,fe  denote  the  bivariate  normal  density  function  with  parameters 
(j,  , a , a,  a , and  p (sometimes  denoted  X,,  X?,  X.,  X.  , X-  for  con- 
venience)  by  0(x,y).  Then,  in  Case  I,  the  likelihood  of  a sample  of  n 
independent  observations  all  in  a region  R is 


4 

p i=l 


where 


p = Pr/”(x,y)  in  R J = 0(x,y)  dx  dy  j 

A 

and,  in  Case  II,  the  likelihood  of  a sample  of  N independent  observa 
tions  of  which  n observations  occur  in  R and  N-n  elsewhere  is 
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The  partial  derivative  of  the  logarithmic  likelihood, L , with  respect  to 
one  of  the  parameters,  say  X,  is 
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3x 


l0e 


where  f(n,p)  * n/p  in  Case  I and  f(n,p)  * (N-n)/(l-p)  in  Case  II.  To 
obtain  the  maximum  likelihood  estimates,  all  five  partial  derivatives 
are  equated  to  zero  and  solved  for  the  unknown  parameters. 

Iterative  solution  of  the  maximum  likelihood  equations.  To  solve 
the  five  estimating  equations  simultaneously,  we  propose  a Newton  itera- 
tive procedure.  Choosing  an  initial  trial  solution,  we  approximate  the 
system  of  equations  by  a linear  system  using  the  linear  terms  in  a Taylor 
series  expansion.  Thus 
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where  a subscript  (1)  denotes  evaluation  at  the  first  trial  point  and  a 
superscript  (1)  denotes  the  first  trial  value.  In  matrix  notation  we  have 
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where  / is  the  (column)  vector  with  elements  (;}*1, . . . ,5),  d ^ is 

the  (column)  vector  with  elements  (A.-\.^),  and 
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Second  trial  values  are  obtained  from  the  first  from 


\ij1  lw 


(assuming  non- singular),  and  by  substituting  these  values  for  the 
initial  ones  further  estimates  are  obtained,  and  so  on  until  stability 
is  reached.  (It  may  not  be  necessary  to  recalculate  the  A matrix  at 
each  step  if  its  elements  are  sufficiently  stationary.  T,hen  it  is  recal- 
culated, its  inverse  may  be  obtained  quickly  by  iteration  £~1  J.) 


decision  of  Estimates.  The  asymototic  covariance  matrix  is  the 
inverse  of  the  matrix  with  elements  , and  thus  may  be  estimated 

by  A~\  This  estimate  is  obtained  simultaneously  with  the  solution  of  the 


estimating  equations. 
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Rectangular  truncation.  Here  we  develop  explicitly  the  estimating 
equations  for  a rectangularly  truncated  population.  Let  the  region  R be  a 
rectangle,  bounded  by  the  lines  x ■ hp  x ■ hg,  y * kp  y « kg  (h^<hg,  k^<  kg). 
Then 

k2  h2 

P a \ ^ 0(x,y)  dx  dy  . 
kl  hl 

Wow  0(x,y)  may  be  expressed  as  a power  series  in  p with  Hermite 
functions  as  coefficients 
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where 


0 (t) 
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1 dv 
/?n  dtv 


( v=l,2, . . , ) 
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(Hee  R.  A.  Fisher's  introduction,  pp.xxvi-xxviii,  £" 2j. ) For  negative  sub 
scripts,  the  Hermite  functions  are  defined  by  the  two  following  relations? 
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n .(t)  - 1 - \ — e'*^Jdx 
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t G (t)  - v G ,(t)  + G ..(t)  . 

V v-1  ' v+1 


Since  the  series  in  p converges  uniformly  in  x and  y,  and  since  ^ Gv(t)  dt 

- G ,(t),  we  have 
v-1 


where 
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After  calculating  the  derivatives  of  the  2 function,  we  obtain  the 
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following  derivatives- 
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all  others  being  obtained  by  symmetry,  ('i'he  interchange  of  x and  y re- 
quires the  interchange  of  the  order  of  the  subscripts  on  the  2 functions.) 

r,s 

In  Case  I, 


(3) 


f(n,p)  = - n/p2 


and  in  Case  IT, 


(U) 


gp  f(n,p)  * (It-n)/(l-p)2  . 


Calculating  the  derivatives  of  log  0(x,y),  we  find' 
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all  others  being  obtained  by  symmetry;  we  have  denoted 
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i 2 (51?*)  (S®1  (r,s  - 0,1,2)  . 


rs  n . ' a ' ' a 

i=*l  x y 


Using  (1)  and  (5),  we  may  now  calculate  the  elements  of  the  (f  vector: 


It  " mf(n>P)  ^ + Z |x*  102  5 

i*l 


and  using  (1),  (2),  (3),  (U),  and  (6),  we  may  calculate  the  elements  of  the 


A matrix 
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(The  required  2 functions  may  be  computed  from  tables  of  the  tetrachoric 
r jS 

functions  /~3  7,  using  the  relation  /Tt  Tv(t)  = Gv_1(t).  ) 

The  estimates  for  truncation  over  an  infinite  quarter-plane  may  be 
obtained  by  letting  one  of  the  h‘s  and  one  of  the  k’s  in  the  discussion 
above  go  to  + co. 

Initial  estimates  may  be  obtained  by  approximating  the  region  R by 
an  infinite  strip  and  using  the  linear  truncation  estimation  method  which 
follows . 


Linear  truncation.  We  shall  consider  independently  the  estimation 
problem  when  the  region  R is  an  infinite  strip,  bounded  by  x « h]_  and  x * h2 
(h^<  h^).  Here,  we  shall  distinguish  three  ones*  (I)  the  number  of  un- 
measured observations  is  unknown,  (II)  the  numbers  of  unmeasured  observations 
in  each  truncated  half -plane  are  known  (say  n^  observations  in  R^  = 

/”(x,y)’  x<h^7  and  n2  in  Rg  = /“(x,y)  x>h2_7,  n-,  + n?  * N-n),  and  (III) 
only  the  total  number  N-n  of  unmeasured  observations  is  known. 

Now  the  marginal  distribution  of  x is  independent  of  a , cr  , and  p, 

V J 

and  is  siimoly  a truncated  univariate  normal  distribution.  Thus,  in  all 
three  cases,  |ix  and  a ■ may  be  estimated  by  the  methods  of  A.  C.  Cohen  £hj 
for  truncated  univariate  normal  distributions,  the  three  cases  above  corres- 
ponding to  the  three  cases  enumerated  by  Cohen. 

The  likelihood  functions  in  the  three  cases  are  ■ 
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where 

p * Pr/~ (x,y)  In  R_7  - K^)  - $<£j£*) 

./V  A 

Pi  = Pr/~(x,y)  in  R^.7  * £(^5^) 

P2  ■ pr/“(x,y)  in  R ZJ  - 1 - 


t 

i (t)  - ( -i-  e“t2//2  dt  ; 
) /2n 
-00 


l<i  and  are  constants.  Since  p,  p^  and  are  aIl  independent  of  |xy,  ay, 
and  p , the  maximum  likelihood  equations  for  these  three  parameters  will  be 
the  same  in  all  three  cases:  namely,  from  (5)  and  (9), 
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where  L denotes  the  logarithmic  likelihood,  m being  defined  by  (7).  Hav- 
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ing  obtained  estimates  of  | ± and  cr  by  Cohen's  method,  these  three  equations 

X X 

may  be  solved  for  estimates  of  (i  , a , and  p , yielding  ( after  some  algebraic 

y y 

manipulation) : 
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where  m ' denotes  sample  moments  about  the  origin.  If  u * m,  then  we 
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substitute 
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in  equations  (10)  and  (11). 

Since  the  estimates  of  p,  and  a are  independent  of  p , o , and  p, 

xx  y y 

we  find  that  the  A matrix  is  now  the  direct  sum  of  two  sub-matrices,  and 
likewise  for  its  inverse.  The  first  inverse  sub-matrix  (corresponding  to 
p and  a ) is  given  by  Cohen  /~U  7>  the  second  may  be  obtained  by  inverting 
the  matrix  whose  elements  are  given  by  (8). 
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